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Strongly correlated Fermi systems are among the most intriguing, best experimentally studied 
and fundamental systems in physics. There is, however, lack of theoretical understanding in this 
field of physics. The ideas based on the concepts like Kondo lattice and involving quantum and 
thermal fluctuations at a quantum critical point have been used to explain the unusual physics. 
Alas, being suggested to describe one property, these approaches fail to explain the others. This 
means a real crisis in theory suggesting that there is a hidden fundamental law of nature. It turns 
out that the hidden fundamental law is well forgotten old one directly related to the Landau — Migdal 
quasiparticles, while the basic properties and the scaling behavior of the strongly correlated systems 
can be described within the framework of the fermion condensation quantum phase transition 
(FCQPT). The phase transition comprises the extended quasiparticle paradigm that allows us to 
explain the non-Fermi liquid (NFL) behavior observed in these systems. In contrast to the Landau 
paradigm stating that the quasiparticle effective mass is a constant, the effective mass of new 
quasiparticles strongly depends on temperature, magnetic field, pressure, and other parameters. Our 
observations are in good agreement with experimental facts and show that FCQPT is responsible for 
the observed NFL behavior and quasiparticles survive both high temperatures and high magnetic 
fields. 
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1. INTRODUCTION 

A. B. Migdal's contribution to modern theoretical 
physics is very impressive. His endowment in and deep 
understanding of different domains of physics, including 
the nuclear and many-body physics, based on Fermi- 
liquid approach, is outstanding. In Migdal seminal 
papers a solid base for studying strongly interacting 
Fermi systems and phase transitions occurring in them 
has been established [l|, |2|. Migdal's daring ideas of 
phase transitions related to tt condensation in nuclei 
and neutron stars Q inspired a theory of fermion 
condensation that has permitted to construct a new 
class of Fermi liquids, new quasiparticles and new type 
of merging of single-particle levels of both finite and 
infinite Fermi systems like nuclear, atomic and solid 
state systems ','■'> 7|. The new class of Fermi liquids is 
represented by strongly correlated Fermi systems where 
enormous number of experimental facts are collected. 
Understanding the physics of these systems stimulates 
intensive studies of the possible manifestation of fermion 
condensation in other areas, as it has happened in 
the case of metal superconductivity, whose ideas were 
successfully used in describing atomic nuclei [l| and 
in a possible explanation of the origin of the mass of 
elementary particles. Therefore, we expect that the ideas 
associated with the new fermion condensation quantum 
phase transition [8| in one area of research stimulates 
intensive studies of the possible manifestation of such a 
transition in other areas. 
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Strongly correlated Fermi systems represented by 
heavy fermion (HF) metals and quasi-two-dimensional 
3 He are among the most intriguing, best experimentally 
studied and fundamental systems in physics, which until 
very recently have lacked theoretical explanations @. 
These are also a field never far from applications in 
synthesis of novel materials for cryogenics, rare earth 
magnets and applied superconductivity. The properties 
of these materials differ dramatically from those of 
ordinary Fermi systems [cj-[l4]|. Their behavior is so 
unusual that the traditional Landau quasiparticles 
paradigm does not apply to it. The paradigm states 
that the properties is determined by quasiparticles 
whose dispersion is characterized by the effective mass 
M* which is independent of temperature T, the 
number density x, magnetic field B and other external 
parameters. The above systems are, however, in defiance 
of theoretical understanding. The ideas based on the 
concepts (like Kondo lattice involving quantum and 
thermal fluctuations at a quantum critical point (QCP) 
have been used to explain the unusual physics of these 
systems known as non-Fermi liquid (NFL) behavior. 
Alas, being suggested to describe one property these 
approaches fail to explain the others. This means 
a real crisis in theory suggesting that there is a 
hidden fundamental law of nature, which remains to 
be recognized. It is widely believed that utterly new 
concepts are required to describe the underlying physics. 
There is a fundamental question: how many concepts do 
we need to describe the above physical mechanisms? This 
cannot be answered on purely experimental or theoretical 
grounds. Rather, we have to use both of them. For 
instance, in the case of metals with heavy fermions, the 
strong correlation of electrons leads to a renormalization 
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of the effective mass of quasiparticles, which may exceed 
the ordinary, "bare mass by several orders of magnitude 
or even become infinitely large at temperatures T — > 0. 
Moreover, the effective mass strongly depends on the 
temperature, pressure, or applied magnetic field. Such 
metals exhibit NFL behavior and unusual power laws 
of the temperature dependence of the thermodynamic 
properties at low temperatures. 

The Landau theory of the Fermi liquid has remarkable 
results in describing a multitude of properties of the 
electron liquid in ordinary metals. Fermi liquids of the 
3 He type and nuclear liquid |15l - [l8| . The theory is based 
on the assumption that elementary excitations determine 
the physics at low temperatures. These excitations 
behave as quasiparticles, have a certain effective mass, 
and, judging by their basic properties, belong to the class 
of quasiparticles of a weakly interacting Fermi gas. Hence, 
the effective mass M* is independent of the temperature, 
pressure, and magnetic field strength and is a parameter 
of the theory. The Landau Fermi liquid (LFL) theory 
fails to explain the results of experimental observations 
related to the dependence of M* on the temperature 
T, magnetic field B, pressure, etc.; this has led to the 
conclusion that quasiparticles do not survive in strongly 
correlated Fermi systems and that the heavy electron 
does not retain its identity as a quasiparticle excitation, 
see e.g. 0, EI HI- 

The unusual properties and NFL behavior observed 
in high-T c superconductors, HF metals and 2D Fermi 
systems are assumed to be determined by various 
magnetic quantum phase transitions [91413. EH, [2l| . 
Indeed, when reasoning by analogy with respect to the 
second order phase transitions, one can assume that 
a phase transition responsible for the NFL behavior 
taking place up to lowest accessible temperatures is 
located at T = 0. Since a quantum phase transition 
occurs at T = 0, the control parameters are the 
composition, electron (hole) number density x, pressure, 
magnetic field strength B, etc. A quantum phase 
transition occurs at a quantum critical point, which 
separates the ordered phase that emerges as a result of 
quantum phase transition from the disordered phase. It 
is usually assumed that magnetic (e.g., ferromagnetic 
and antiferromagnetic) quantum phase transitions are 
responsible for the NFL behavior. The critical point of 
such a phase transition can be shifted to absolute zero 
by varying the above parameters. 

Universal behavior can be expected only if the system 
under consideration is very close to a quantum critical 
point, e.g., when the correlation length is much longer 
than the microscopic length scale, and critical quantum 
and thermal fluctuations determine the anomalous 
contribution to the thermodynamic functions of strongly 
correlated Fermi system. Quantum phase transitions of 
this type are so widespread [IM1 Q El 

that we call 

them ordinary quantum phase transitions [22l |. In this 
case, the physics of the phenomenon is determined by 
thermal and quantum fluctuations of the critical state, 



while quasiparticle excitations are destroyed by these 
fluctuations. Conventional arguments that quasiparticles 
in strongly correlated Fermi liquids "get heavy and die "at 
a quantum critical point commonly employ the well- 
known formula based on the assumptions that the z- 
factor (the quasiparticle weight in the single-particle 
state) vanishes at the points of second-order phase 
transitions [2fJ. However, it has been shown that this 
scenario is problematic (23l . [24| . 

The fluctuations in the order parameter developing 
an infinite correlation length and the absence of 
quasiparticle excitations are considered as the main 
reason for the NFL behavior of heavy- fermion metals, 2D 
fermion systems and high-T c superconductors [ill EH, EH 
EslHHl- This approach faces certain difficulties, however. 
Critical behavior in experiments with metals containing 
heavy fermions is observed at high temperatures 
comparable to the effective Fermi temperature Tfc. 
For instance, the thermal expansion coefficient a(T), 
which is a linear function of temperature for normal 
LFL, a(T) oc T, demonstrates the \JT temperature 
dependence in measurements involving CeNi2Ge2 as the 
temperature varies by two orders of magnitude (as it 
decreases from 6 K to at least 50 mK) [21]. Such 
behavior can hardly be explained within the framework 
of the critical point fluctuation theory. Obviously, such 
a situation is possible only as T — > 0, when the critical 
fluctuations make the leading contribution to the entropy 
and when the correlation length is much longer than 
the microscopic length scale. At a certain temperature 
Tfc, this macroscopically large correlation length must 
be destroyed by ordinary thermal fluctuations and the 
corresponding universal behavior must disappear. 

In the rest of this paper, we show that the fermion 
condensation quantum phase transition (FCQPT) Q 
is indeed responsible for the observed fascinating NFL 
behavior of strongly correlated Fermi systems and 
quasiparticles survive both high temperatures and high 
magnetic fields. In Section [2] we give a detailed 
consideration of experimental evidences in favor of 
existence of quasiparticles, and formulate both a scaling 
behavior of strongly correlated Fermi systems and the 
extended quasiparticle paradigm. Then in Section [31 
we consider the properties of Landau Fermi liquid. In 
Section we demonstrate that the Landau equation for 
the effective is not a phenomenological one and can 
be derived using the methods of Density Functional 
Theory. Thus, we establish the extended quasiparticle 
paradigm. FCQPT and a phase diagram of heavy fermion 
system located in the vicinity of FCQPT are investigated 
in Section [5] We propose that the phase diagram of 
systems located near FCQPT are strongly influenced by 
control parameters such as a chemical pressure, pressure 
or magnetic field. We find that under the application 
of the chemical pressure (positive/negative) QCP is 
destroyed or converted into a quantum critical line, 
correspondingly. In Section [S] we establish that heavy 
fermion quasiparticles do exist in a very wide range of 
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both temperatures T and magnetic fields B. Finally, in 
Section [7] our results are summarized and discussed. 



2. EXPERIMENTAL HINTS AT SCALING 
BEHAVIOR AND QUASIPARTICLES 

There is a difficulty in explaining the restoration of the 
LFL behavior under the application of magnetic field B, 
as observed in HF metals and in high-T c superconductors 
d, HE H3- For the LFL state as T -> 0, the electric 
resistivity p(T) = po + AT 2 , the heat capacity C(T) = 
7oT, and the magnetic susceptibility \ = const. It turns 
out that the coefficient A(B), the Sommerfeld coefficient 
7o(-B) = C/T cx M*, and the magnetic susceptibility 
x(B) depend on the magnetic field strength B such that 
A(B) cx 7o(-B) and A(B) cx x 2 (B), which implies that 
the Kadowaki- Woods relation K = A(B)/~f 2 (B) [1| 
is _B-independent and is preserved |27| . Such universal 
behavior, quite natural when quasiparticles with the 
effective mass M* playing the main role, can hardly 
be explained within the framework of approach that 
presupposes the absence of quasiparticles, which is 
characteristic of ordinary quantum phase transitions in 
the vicinity of QCP. Indeed, there is no reason to expect 
that 7o, x and A are affected by the fluctuations in a 
correlated fashion. 

For instance, the Kadowaki— Woods relation does not 
agree with the spin density wave scenario [27i | and 
with the results of research in quantum criticality based 
on the renormalization-group approach [29]. Moreover, 
measurements of charge and heat transfer have shown 
that the Wiedemann— Franz law holds in some high-T c 
superconductors [H,|30| and HF metals [311-I3411 . All this 
suggests that quasiparticles do exist in such metals, and 
this conclusion is also corroborated by photoemission 
spectroscopy results [35l l36j. 

The inability to explain the behavior of heavy-fermion 
metals while staying within the framework of theories 
based on ordinary quantum phase transitions implies 
that another important concept introduced by Landau, 
the order parameter, also ceases to operate. Thus, we are 
left without the most fundamental principles of many- 
body quantum physics (TH HH, [l8| , and many interesting 
phenomena associated with the NFL behavior of strongly 
correlated Fermi systems remain unexplained. 

The NFL behavior manifests itself in the power- 
law behavior of the physical quantities of strongly 
correlated Fermi systems located close to their QCPs, 
with ex pon ents different from those of a Fermi liquid [9- 
0, [H, IMI, 137^39). It is common belief that the main 
output of theory is the explanation of these exponents 
which are at least depended on the magnetic character 
of QCP and dimensionality of the system. On the 
other hand, the observed behavior of the thermodynamic 
properties cannot be captured by these exponents as seen 
from Figs. [T]— [31 The behavior of the entropy S(T) of two- 
dimensional (2D) 3 He [40(| shown in Fig. Q] is positively 
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1: The entropy of 2D 3 He as a function of the density x 
versus T [4p| |. The density is depicted in the legend. 



different from that described by a simple function a\T a ' 1 
where a\ is a constant and is the exponent. It is seen 
from Fig. Q] that at the low densities x ~ 7 nm -2 the 
entropy demonstrates the LFL behavior characterized 
by a linear function of T with a% = 1. The behavior 
becomes quite different at the higher densities at which 
S(T) has an inflection point. Obviously, at the inflection 
point S(T) cannot be fit by the simple function aiT a2 . 
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Phc. 2: Electronic specific heat of YbRh2Si2, C/T, versus 
temperature T as a function of magnetic field B [37} shown 
in the legend. 

As seen from Fig.[2j the specific heat C/T measured on 
YbRli2Si2 exhibits a behavior that is to be described 
function of both temperature T and magnetic B 
field rather than by a single exponent. One can see that 
at low temperatures C/T demonstrates the LFL behavior 
which is changed by the transition regime at which C /T 
reaches its maximum and finally C/T decays into NFL 
behavior as a function of T at fixed B. It is clearly seen 
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from Fig. [5] that, in particularly in the transition regime, 
these exponents may have little physical significance. 
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Phc. 3: The normalized effective mass = M* /M^.j versus 
normalized temperature Tjv = T/Tm- is extracted from 
the measurements of the specific heat C/T on YbRh2Si2 [Hj] 
shown in Fig. [3] The values of the field B is listed in the 
legend. 

Figure [3] displays measurements of C/T on YbRh2Si2 
[38| in wide range of both temperature T and magnetic 
field B variations, ft is seen that both the NFL behavior 
and the range extends at least up to twenty Kelvins and 
up to 18 T. Figure [3] demonstrates that the LFL state at 
low temperatures and NFL one at higher temperatures 
are separated by the transition regime at which C /T 
reaches its maximum value M M (B). Thus, we again 
conclude that the observed NFL behavior exhibiting 
the maximum and extending to high temperatures can 
be hardly explained in the framework of the ordinary 
quantum phase transitions interpreting the exponents. 

In order to show that the behavior of S and C/T 
displayed in Figs. HH3] is of generic character, we 
remember that in the vicinity of QCP it is helpful 
to use "internal "scales to measure the effective mass 
M* (x C/T and temperature T" 0, EE S2|- The internal 
scales of the thermodynamic functions such as S or 
C/T are related to "peculiar points "like the inflection 
or maximum. Since the entropy has no maxima, its 
normalization is to be performed in the inflection point 
taking place at T = Ti n f. Note that Tin/ is a function 
of x, it is seen from Fig. [T] that the inflection point 
moves towards lower temperatures at elevating x. The 
normalized entropy Sn as a function of the normalized 
temperature TV = T/Ti n t = y is reported in Fig. [4] 
We normalize the entropy by its value at the inflection 
point S N (y) = S(y)/S(l). As seen from Fig. gj the 
normalization revels the scaling behavior of Sn, that 
is the curves at different temperatures and densities 
x merge into a single one in terms of the variable y. 
We have excluded the experimental data taken at x < 
8 nm -2 since the corresponding curves do not contain 
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Phc. 4: The normalized entropy Sn, extracted from the 
experimental facts displayed in Fig. U as a function of x 
shown in the legend versus normalized Tat. The arrow shows 
the inflection point. 



the inflection points. It is seen from Fig. U that <Sjv(y) 
extracted from the measurements is not a linear function 
of y, as would be for a LFL, and shows the scaling 
behavior over three decades in normalized temperature 
V- 

As seen from Fig. [51 a maximum structure in C /T oc 
M M at temperature Tm appears under the application 
of magnetic field B and Tm shifts to higher T as B 
is increased. The value of the Sommerfeld coefficient 
C/T = 70 is saturated towards lower temperatures 
decreasing at elevated magnetic field. To obtain the 
normalized effective mass M N , we use M M and Tm as 
"internal"scales: The maximum structure in C/T was 
used to normalize C/T, and T was normalized by Tm- 
In Fig. [5] the normalized effective mass — M*/M M 
as a function of normalized temperature Tn — T/Tm is 
shown by geometrical figures. Note that we have excluded 
the experimental data taken in magnetic field B = 0.06 
T. In that case, Tm — > and the corresponding Tm and 
M M are unavailable. It is seen that the LFL state and 
NFL one are separated by the transition regime at which 
M N reaches its maximum value. Figure [5] reveals the 
scaling behavior of the normalized experimental curves: 
The curves at different magnetic fields B merge into a 
single one in terms of the normalized variable y = T/Tm- 
As seen from Fig.[5l the normalized effective mass M N (y) 
extracted from the measurements is not a constant, as 
would be for a LFL, and shows the scaling behavior 
over three decades in normalized temperature y. It is 
seen from Figs. [2] and [5] that the NFL behavior and the 
associated scaling extend at least to temperatures up to 
few Kelvins. 

In order to get a deep insight in the behavior of C/T 
displayed in Fig. [3J we again use the internal scales as 
it was done when constructing M N displayed in Fig. [5] 
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Phc. 5: 

temperature Tjv. is extracted from the measurements of 
the specific heat C/T on YbRli2Si2 in magnetic fields B (37J 
listed in the legend. Constant effective mass M£ inherent in 
normal Landau Fermi liquids is depicted by the solid line. 



As a result, Fig. [5] reveals the scaling behavior of the 
normalized experimental curves: The curves at different 
magnetic fields B merge into a single one in terms of the 
normalized variable y = T/Tm- As seen from Fig. [6j the 
normalized effective mass M N (y) is not a constant, as 
would be for a LFL, and shows the scaling behavior over 
two decades in normalized temperature y. It is seen from 
Figs. [3] and [6] that the NFL behavior and the associated 
scaling extend at least both to high temperatures up to 
twenty Kelvins and magnetic fields up to 18 T. 
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Phc. 6: The normalized effective mass — M* /M^ versus 
normalized temperature Tjv = T/Tm. is extracted from 
the measurements of the specific heat C/T on YbRli2Si2 HH 
shown in Fig. The values of the field B is listed in the 
legend. 



Taking into account the scaling behavior revealed 
in Figs. we conclude that a challenging problem 



for theories considering the NFL behavior of strongly 
correlated Fermi liquids is to explain both the scaling and 
the shape of the thermodynamic functions like Sn an d 
M N . While the theories calculating only the exponents 
that characterize Sn and M N at y 3> 1 deal with a part 
of the observed facts related to the problem and overlook, 
for example, consideration of the transition and LFL 
regimes. Another part of the problem is the remarkably 
large ranges of temperature, magnetic field and the 
density over which the NFL behavior and the scaling 
are observed. Concepts based on the Kondo lattice and 
scenarios where fluctuations in the order parameter is 
sufficiently big and the correlation time is sufficiently 
large to develop the NFL behavior can hardly match up 
such high temperatures and explain the observed scaling 
behavior. 

As we will see below, the large temperature ranges are 
precursors of new quasiparticles, and it is the scaling 
behavior of the normalized effective mass that allows 
us to explain the thermodynamic of HF metals at the 
transition and NFL regimes. Taking into account the 
simple behavior shown in Figs. UJ— [51 we ask the question: 
what theoretical concepts can replace the Fermi-liquid 
paradigm with the notion of the effective mass in cases 
where Fermi-liquid theory breaks down? To date such 
a concept is not available [111]. Therefore, we focus 
on the concept of FCQPT preserving quasiparticles 
and intimately related to the unlimited growth of 
M*. It was shown that it is capable of revealing the 
scaling behavior of the effective mass and delivering 
an adequate theoretical explanation of a vast majority 
of experimental results in different HF metals @. 
In contrast to the Landau paradigm based on the 
assumption that M* is a constant as shown by the 
solid line in Fig. [5J in FCQPT approach the effective 
mass M* of new quasiparticles strongly depends on 
T, x, B etc. Therefore, in accord with numerous 
experimental facts the extended quasiparticles paradigm 
is to be introduced. The main point here is that 
the well-defined quasiparticles determine as before the 
thermodynamic, relaxation and transport properties of 
strongly correlated Fermi-systems in large temperature 
ranges, while M* becomes a function of T. x, B etc. 
The FCQPT approach had been already successfully 
applied to describe the thermodynamic properties of 
such different strongly correlated systems as 3 He on the 
one hand and complicated HF compounds on the other 

H Hi HE 51. 



Since we are concentrated on properties that are non- 
sensitive to the detailed structure of the system we avoid 
difficulties associated with the anisotropy generated 
by the crystal lattice of solids, its special features, 
defects, etc., We study the universal behavior of strongly 
correlated Fermi-systems located near their QCP at low 
temperatures using the model of a homogeneous HF 
liquid [U, The model is quite meaningful because 
we consider the scaling behavior exhibited by these 
materials at low temperatures Q . The scaling properties 
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of the normalized effective mass that characterizes them, 
are determined by momentum transfers that are small 
compared to momenta of the order of the reciprocal 
lattice length. The high momentum contributions can 
therefore be ignored by substituting the lattice for the 
jelly model. While the values of the scales like the 
maximum Mj^ of the effective mass and Tm at which 
takes place are determined by a wide range of 
momenta and thus these scales are controlled by the 
specific properties of the system. On the other hand, 
the dependencies of these scales on magnetic fields, 
temperature, the density of system, pressure etc are again 
controlled by small momenta and can be analyzed within 
the model of HF liquid. 



3. NORMAL FERMI LIQUIDS 

One of the most complex problems of modern 
condensed matter physics is the problem of the structure 
and properties of Fermi systems with large inter particle 
coupling constants. Theory of Fermi liquids, later called 
"normal was first proposed by Landau as a means 
for solving such problems by introducing the concept 
of quasiparticles and amplitudes that characterize the 
effective quasiparticle interaction (TBI . [TfjJ. The Landau 
theory can be regarded as an effective low-energy theory 
with the high-energy degrees of freedom eliminated by 
introducing amplitudes that determine the quasiparticle 
interaction instead of the strong inter particle interaction. 
The stability of the ground state of the Landau Fermi 
liquid is determined by the Pomeranchuk stability 
conditions: stability is violated when at least one Landau 
amp litude becomes negative and reaches its critical value 
. We note that the new phase in which stability 
is restored can also be described, in principle, by the LFL 
theory. 

We begin by recalling the main ideas of the LFL 
theory [l8j]. The theory is based on the 

quasiparticle paradigm, which states that quasiparticles 
are elementary weakly excited states of Fermi liquids and 
are therefore specific excitations that determine the low- 
temperature thermodynamic and transport properties of 
Fermi liquids. In the case of the electron liquid, the 
quasiparticles are characterized by the electron quantum 
numbers and the effective mass M*. The ground state 
energy of the system is a functional of the quasiparticle 
occupation numbers (or the quasiparticle distribution 
function) n(p,T), and the same is true of the free 
energy F(n(p,T)), the entropy 5(n(p, T)), and other 
thermodynamic functions. We can find the distribution 
function from the minimum condition for the free energy 
F = E — TS (here and in what follows ks = h = 1) 



Here fi is the chemical potential fixing the number density 

dp 



. = /»(p,T)pg 



and 



<P,T) 



SE(n(p,T)) 
Sn(p,T) 



(2) 



(3) 



is the quasiparticle energy. This energy is a functional of 
n(p,T), in the same way as the energy E is: e(p,T,n). 
The entropy S(n(p,T)) related to quasiparticles is given 
by the well-known expression (IH 

S(n(p,T)) = -2 An(p,r)ln(n(p,T)) + (l-n(p,T)) 

dp 



ln(l-n(p,T))] 



(2tt) 3 ' 



(4) 



which follows from combinatorial reasoning. Equation (fTJ) 
is usually written in the standard form of the Fermi — 
Dirac distribution. 



n(p,T) 



1 + exp 



(e(p,T)-/z) 



T 



(5) 



At T — >• 0, (H} and (5| have the standard solution 
n(p, T —> 0) — s- 9(pf —p) if the derivative de(p ~ pF)/dp is 
finite and positive. Here pp is the Fermi momentum and 
8(pf — p) is the step function. The single particle energy 
can be approximated as e(p ~ p-p) — ~ Pf(p — Pf)/Mj j , 
and M£ inversely proportional to the derivative is the 
effective mass of the Landau quasiparticle, 



1 



1 de(p,T = 0) 
p dp 



(6) 



In turn, the effective mass Mi is related to the bare 



electron mass 

[HEIDI 



Mi 



by the well-known Landau equation 



1 

m 
On 



E 

CTl 



Pf 

CT1 (pi,T) dpi 



PfPi „ , , 

3 — - f V,<Ti(.PF,PlJ 



dpi 



(2tt) 3 



(7) 



where ^^(pFiPx) i s the Landau amplitude, which 
depends on the momenta pp and p and the spins a. For 
simplicity, we ignore the spin dependence of the effective 
mass, because M£ is almost completely spin-independent 
in the case of a homogeneous liquid and weak magnetic 
fields. The Landau amplitude F is given by 



F<j,o l (p,Pi,n) 



6 2 E(n) 



Sn a (p)Sn ai (pi)' 



(8) 



The stability of the ground state of LFL is determined 
by the Pomeranchuk stability conditions: stability is 
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violated when at least one Landau amplitude becomes 
negative and reaches its critical value fla [H, E3| 



■(2L + 1). 



(9) 



Here F£ and F[ are the dimensionless spin-symmetric 
and spin- antisymmetric Landau amplitudes, L is 
the angular momentum related to the corresponding 
Legendre polynomials Pj,, 



F(pa, Pl <r x ) = 1 £ P L (6) CTl + F£] . 



(10) 



L=0 



Here is the angle between momenta p and pi and the 
density of states N — M£pp/(27r 2 ). It follows from Eq. 
© that 



Mr Ff 

_ L = l + -L. 

m 3 



(11) 



In accordance with the Pomeranchuk stability conditions 
it is seen from Eq. (fTTj) that Ff > —3, otherwise the 
effective mass becomes negative leading to unstable state 
when it is energetically favorable to excite quasiparticles 
near the Fermi surface. 



4. EFFECTIVE MASS AND SCALING 
BEHAVIOR 

It is common belief that the equations of the Section [3] 
are phenomenological and inapplicable to describe Fermi 
systems characterized by the effective mass M* strongly 
dependent on temperature, external magnetic fields B, 
pressure P etc. To derive the equation determining the 
effective mass, we consider the model of a homogeneous 
HF liquid and employ the density functional theory 
for superconductors (SCDFT) [45| which allows us to 
consider £ as a functional of the occupations numbers 
n (p) O l46U48j . As a result, the ground state energy 
of the normal state E becomes the functional of the 
occupation numbers and the function of the number 
density x, E = E(n(p),x), while Eq. ([3]) gives the 
single-particle spectrum. Upon differentiating both sides 
of Eq. with respect to p and after some algebra and 
integration by parts, we obtain [23l [ill . |46| . |47| 

de(p) = £ 
dp m 

To calculate the derivative de(p)/dp, we employ the 
functional representation 

/rp dp 
7^-n(p)—-^ 
Ira (27r) J 



f . N 9n(pi) dpi 

]F(p, Pl ,n) ^ (27r)3 - (12) 



It is seen directly from Eq. ([12]) that the effective mass 
is given by the well-known Landau equation 



1 J_ , f PfPi 

M* ~ m 



f PfPi „, ,dn(pi) dpi f . 

/ — — ^(PF,Pi,n)— ^ -^3. (14 

J p% dpx 2tt) 3 



For simplicity, we ignore the spin dependencies. To 
calculate M* as a function of T, we construct the free 
energy F — E — TS, where the entropy S is given by 
Eq. Minimizing F with respect to n(p), we arrive 
at the Fermi — Dirac distribution, Eq. ([5]). Due to the 
above derivation, we conclude that Eqs. (fl"2"T) and ([M)) are 
exact ones and allow us to calculate the behavior of both 
cte(p) / dp and M* which now is a function of temperature 
T, external magnetic field B, number density x and 
pressure P rather than a constant. As we will see it is 
this feature of M* that forms both the scaling and the 
NFL behavior observed in measurements on HF metals. 

In LFL theory it is assumed that M£ is positive, finite 
and constant since the integral on the right hand side 
of Eq. (fT4| represents a small correction in the case of 
normal metals. As a result, the temperature-dependent 
corrections to Af£, the quasiparticle energy e(p) and 
other quantities begin with the term proportional to T 2 
in 3D systems and with the term proportional to T in 2D 
one [49| . The effective mass is given by Eq. ([7]), and the 
specific heat C is [15] 



2^NT 8S 
C=^— =7oT = T— , 



and the spin susceptibility 



X 



^{i + F%y 



(15) 



(16) 



where /is is the Bohr magneton and 70 oc Af£. In the 
case of LFL the electrical resistivity at low T is given by 

m 



P (T) = po + AT 2 



(17) 



where po is the residual resistivity and A is the coefficient 
determining the charge transport. The coefficient A is 
proportional to the quasiparticle-quasiparticle scattering 
cross-section. Equation (|17j) symbolizes and defines the 
LFL behavior observed in normal metals. 

Equation (fLT| at T — 0, combined with the fact that 
n(p,T = 0) becomes 9(pp — p), yields the well-known 
result 150-521 



M* 
m 



1 



1 - /A/3 



where fjj 



Nof 1 , N — mpF/(27r 2 ) is the density 
of states of a free Fermi gas and ^(pf^Pf) is the p- 
wave component of the Landau interaction amplitude. 
Because x = p|/37r 2 in the Landau Fermi-liquid theory, 
the Landau interaction amplitude can be written as 
/at(PF)Pf) — fh( x )- Provided that at a certain critical 
point XfGj the denominator (1 — fh(x)/3) tends to zero, 
i.e., (1 — fl r (x)/3) oc (x — xfc) + a{x — xfg) 2 + ••• -> 0, 
we find that 153, 541 



M*{x) 



01 



a 2 1 
oc -. 

- xfc r 



(18) 



where a\ and ai are constants and r = (x — xfcV^fc 
is the "distance" from QCP x FC at which M*(x -» 
^fc) — > oo. We note that the divergence of the effective 
mass given by Eq. (|18[) does preserve the Pomeranchuk 
stability conditions for /jy is positive, see Eq. ©. 
Equations (|lip and (|18p seem to be different but it is 



not the case since fl. 



oc m, while oc M* and Eq. (|TT)I 
represents an implicit formula for the effective mass. 

The behavior of M* (x) described by formula ((TBI is in 
good agreement with the results of experiments [55l457j 
and calculations [58-60]. In the case of electron systems, 
Eq. (TTgj) holds for x > xfc, while for 2D 3 He we have 
x < xfc so that always r > [3, HH . Such behavior of 
the effective mass is observed in HF metals, which have 
a fairly flat and narrow conductivity band corresponding 
to a large effective mass, with a strong correlation and 
the effective Fermi temperature Tj~ ~ p F /M*(x) of the 
order of several dozen degrees kelvin or even lower (e.g., 
see Ref. [jj). 

The effective mass as a function of the electron density 
x in a silicon MOSFET (Metal Oxide Semiconductor 
Field Effect Transistor), approximated by Eq. (p~8|) . is 
shown in Fig. [JJ The parameters ai, a 2 and xfc are 
taken as fitting. We see that Eq. fTB"]) provides a good 
description of the experimental results. The divergence 




x [nm" ] 



Phc. 8: The ratio M* /M in 2D 3 He as a function of the density 
x of the liquid, obtained from heat capacity and magnetization 
measurements. The experimental data are marked by squares 
[HL 113 1 an d the solid line represents the function given by 
Eq. 011). 
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Phc 7: The ratio M*/M in a silicon MOSFET as a function 
of the electron number density x. The squares mark the 
experimental data on the Shubnikov-de Haas oscillations, and 
the data obtained by applying a parallel magnetic field are 
marked by circles [561 1571 163 |. The solid line represents the 
function ((18)) . 

of the effective mass M*(x) discovered in measurements 
involving 2D 3 He [HE El] is illustrated in Fig. |SJ 
Figures [JJ and [SJ show that the description provided by 
Eq. (|T8j) does not depend on elementary Fermi particles 
constituting the system and is in good agreement with 
the experimental data. It is instructive to briefly explore 
the scaling behavior of M* in order to illustrate the 



ability of the quasiparticle extended paradigm to capture 
the scaling behavior. Let us write the quasiparticle 
distribution function as ni(p) = n(p,T) — n(p), with 
n(p) being the step function, and Eq. ([T4|) then becomes 



1 



M*(T) 



1 



+ 



PfPi 



-F(Pf,Pi 



dn x {pi,T) dpi 
dpi 



(2tt) 3 
(19) 

At QCP x — > xfc, the effective mass M*(x) diverges 
and Eq. (1191) becomes homogeneous determining M* as 
a function of temperature while the system exhibits the 
NFL behavior. If the system is located before QCP, M* is 
finite, at low temperatures the integral on the right hand 
side of Eq. ([T9")) represents a small correction to 1 /M* and 
the system demonstrates the LFL behavior seen in Figs. 
[5] and [SJ The LFL behavior assumes that the effective 
mass is independent of temperature, M*(T) ~ const, as 
shown by the horizontal line in Fig. [SJ Obviously, the 
LFL behavior takes place only if the second term on 
the right hand side of Eq. (I19[) is small in comparison 
with the first one. Then, as temperature rises the system 
enters the transition regime: M* grows, reaching its 
maximum Mj^ at T = Tm, with subsequent diminishing. 
As seen from Fig. [SJ near temperatures T > Tm the last 
"traces"of LFL regime disappear, the second term starts 
to dominate, and again Eq. (|19p becomes homogeneous, 
and the NFL behavior is restored, manifesting itself in 
decreasing M* as a function of T. 



9 



FERMION CONDENSATION QUANTUM 
PHASE TRANSITION 



is characterized by the invariant [!, [|| [68[ 



As shown in Section |4j the Pomeranchuk stability 
conditions do not encompass all possible types of 
instabilities and that at least one related to the 
divergence of the effective mass given by Eq. (fTSj) 
was overlooked [3j. This type of instability corresponds 
to a situation where the effective mass, the most 
important characteristic of quasiparticles, can become 
infinitely large. As a result, the quasiparticle kinetic 
energy is infinitely small near the Fermi surface and 
the quasiparticle distribution function n(p) minimizing 
E(n(p)) is determined by the potential energy. This leads 
to the formation of a new class of strongly correlated 
Fermi liquids with FC 0K separated from the 

normal Fermi liquid by FCQPT~|22|, |65H67j|. 

It follows from (fT8l) that at T = and as r — > 
the effective mass diverges, M*(r) — > oo. Beyond the 
critical point xpc, the distance r becomes negative and, 
correspondingly, so does the effective mass. To avoid an 
unstable and physically meaningless state with a negative 
effective mass, the system must undergo a quantum phase 
transition at the critical point x = £fc, which, as we 
will see shortly, is FCQPT [H HIH67|. Because the 
kinetic energy of quasiparticles that are near the Fermi 
surface is proportional to the inverse effective mass, the 
potential energy of the quasiparticles near the Fermi 
surface determines the ground-state energy as x — > xpc- 
Hence, at T — a phase transition reduces the energy of 
the system and transforms the quasiparticle distribution 
function. Beyond QCP x = iefc, the quasiparticle 
distribution is determined by the ordinary equation for a 
minimum of the energy functional Q : 



SE(n(p)) 
Sn(p) 



e(p) = m Pi <P< Pf- 



(20) 



Equation (fSUj) yields the quasiparticle distribution 
function tiq(p) that minimizes the ground-state energy 
E. This function found from Eq. (I2TJ1) differs from the 
step function in the interval from pi to pp, where 
< no(p) < 1, and coincides with the step function 
outside this interval. In fact, Eq. (|20p coincides with 
Eq. ^ provided that the Fermi surface at p = p-p 
transforms into the Fermi volume at pi < p < pf 
suggesting that the single-particle spectrum is absolutely 
"flat" within this interval. A possible solution rio(p) of 
Eq. (fSUj) and the corresponding single-particle spectrum 
e(p) are depicted in Fig.O Quasiparticles with momenta 
within the interval (pf—pi) have the same single-particle 
energies equal to the chemical potential /j, and form FC, 
while the distribution nnCp) describes the new state of the 
Fermi liquid with FC [3|-|5j. In contrast to the Landau, 
marginal, or Luttinger Fermi liquids jTo|, which exhibit 
the same topological structure of the Green's function, 
in systems with FC, where the Fermi surface spreads 
into a strip, the Green's function belongs to a different 
topological class. The topological class of the Fermi liquid 
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e(p) = ju 
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Phc. 9: The single-particle spectrum e(p) and the 
quasiparticle distribution function no(p). Because no(p) is a 
solution of Eq. (|20p at T = 0, we have no(p < pi) — 1, 
< no(pi < p < pf) < 1, and no(p > Pf) = 0, while 
e(pi < p < pf) = fj, and the band becomes completely flat. 
The Fermi momentum pp satisfies the condition pi < pf < Pf- 



N = tr 



f —-G(iu}, p)diG~ 1 (iu}, p), (21) 
j c 



where "tr" denotes the trace over the spin indices of 
the Green's function and the integral is taken along 
an arbitrary contour C encircling the singularity of the 
Green's function. The invariant iV in (f2"Tj) takes integer 
values even when the singularity is not of the pole type, 
cannot vary continuously, and is conserved in a transition 
from the Landau Fermi liquid to marginal liquids and 
under small perturbations of the Green's function. As 
shown by Volovik [ESUllj the situation is quite different 
for systems with FC, where the invariant N becomes 
a half-integer and the system with FC transforms into 
an entirely new class of Fermi liquids with its own 
topological structure. 



Phase diagram of Fermi system with FCQPT 

We start with visualizing the main properties of 
FCQPT. To this end, again consider SCDFT. SCDFT 
states that the thermodynamic potential $ is a universal 
functional of the number density n(r) and the anomalous 
density (or the order parameter) /t(r, ri), providing 
a variational principle to determine the densities. 
At the superconducting transition temperature T c a 
superconducting state undergoes the second order phase 
transition. Our goal now is to construct a quantum phase 
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transition which evolves from the superconducting one, 
since, as seen from Fig. [5J the superconducting order 
parameter = yj n(p)(l — n(p)) is hnite over the 

region (p f ~Pi). 

Let us assume that the coupling constant Ao of 
the BCS-like pairing interaction [69( vanishes, with 
Ao — > making vanish the superconducting gap at 
any finite temperature. In that case, T c — > and the 
superconducting state takes place at T — while at finite 
temperatures there is a normal state. This means that at 
T = the anomalous density 

«(r,r 1 ) = (*t(r)*;(r 1 )) (22) 
is finite, while the superconducting gap 

A(r) = Ao J K{r,r 1 )dv 1 (23) 

is infinitely small [jHIBij . In Eq. ([22]>. the field operator 
^/^(r) annihilates an electron of spin a, a =t, 4- a t t ne 
position r. For the sake of simplicity, we consider the 
model of homogeneous HF liquid 01. Then at T = 0, 
the thermodynamic potential $ reduces to the ground 
state energy E which turns out to be a functional of 
the occupation number n(p) since in that case the order 
parameter n(p) = v(p)u(p) — i/n(p)(l — n(p)). Indeed, 

n(p) = v 2 (p); «(p) = v(p)u(p), (24) 

where u(p) and v(p) are normalized parameters such that 
j^p>) +u 2 (p) — 1 and k(p) = y / n(p)(l — n(p)), see e.g. 

Upon minimizing E with respect to n(p), we obtain 
Eq. (f2"0"|) . As soon as Eq. (|2"0|) has nontrivial solution 
n (p) then instead of the Fermi step, we have < 
no (p) < 1 in certain range of momenta pi < p < pf with 
K(p) = y/no(p)(l — «o(p)) being finite in this range, 
while the single particle spectrum e(p) is flat. Thus, the 
step-like Fermi filling inevitably undergoes restructuring 
and forms FC when Eq. (|20[) possesses for the first time 
the nontrivial solution at x = x c which is QCP of 
FCQPT. In that case, the range vanishes, pi — > Pf — LPf, 
and the effective mass M* diverges at QCP [1, 0, [H| 

M*(x -> x c ) ~ pf dp I™™ ^ U - [2b) 

At any small but finite temperature the anomalous 
density k (or the order parameter) decays and this state 
undergoes the first order phase transition and converts 
into a normal state characterized by the thermodynamic 
potential <&o. Indeed, at T —> 0, the entropy S — 
—d$o/dT of the normal state is given by Eq. (j4|. It is 
seen from Eq. @ that the normal state is characterized 
by the temperature- independent entropy So fl3L l23l Frol | . 
Since the entropy of the superconducting ground state 
is zero, we conclude that the entropy is discontinuous 
at the phase transition point, with its discontinuity 



SS = So- Thus, the system undergoes the first order 
phase transition. The heat q of transition from the 
asymmetrical to the symmetrical phase is q — T c So = 
since T c = 0. Because of the stability condition at 
the point of the first order phase transition, we have 
<&o(n(p)) = 4>(k(p)). Obviously the condition is satisfied 
since q — 0. 

At T — 0, a quantum phase transition is driven by a 
nonthermal control parameter, e.g., the number density 
x. As we have seen in Section [4j at QCP, x — xfg, 
the effective mass diverges. It follows from Eq. (fT5)) that 
beyond QCP, the effective mass becomes negative. To 
avoid an unstable and physically meaningless state with 
a negative effective mass, the system undergoes FCQPT 
leading to the formation of FC. 




Phc. 10: Schematic phase diagram of system with FC. The 
number density x is taken as the control parameter and 
depicted as x/xfc- The dashed line shows M* (x/xfc) as 
the system approaches QCP, x/xfc = 1, of FCQPT which 
is denoted by the arrow. At x/xfc > 1 and sufficiently low 
temperatures, the system is in the LFL state as shown by 
the shadow area. The new location of the system induced 
by the application of the positive pressure is shown by both 
the solid square and the arrow. The vertical arrow illustrates 
the system moving in the LFL-NFL direction along T at 
fixed control parameter. At T — and beyond the critical 
point, x/xfc < 1, the system is at the quantum critical 
line depicted by the dashed line and shown by the vertical 
arrow. The critical line is characterized by the FC state with 
finite superconducting order parameter k. At any finite low 
temperature T > T c = 0, k is destroyed, the system undergoes 
the first order phase transition, possesses finite entropy So and 
exhibits the NFL behavior at any finite temperatures T < Tf . 

A schematic phase diagram of the system which is 
driven to the FC state by variation of x is reported 
in Fig. 1101 Upon approaching the critical density xfc 
the system remains in the LFL region at sufficiently 
low temperatures as it is shown by the shadow area. 
The temperature range of the shadow area shrinks as 
the system approaches QCP, and M*(x/xfg) diverges 
as shown by the dashed line and Eq. (fl"8]l . At QCP xfc 
shown by the arrow in Fig. 1101 the system demonstrates 
the NFL behavior down to the lowest temperatures. 
Beyond the critical point at finite temperatures the 
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behavior remains the NFL and is determined by the 
temperature- independent entropy So (l3l . [70| . In that 
case at T — >• 0, the system is again demonstrating the 
NFL behavior and approaching a quantum critical line 
(QCL) (shown by the vertical arrow and the dashed line 
in Fig. I10p rather than a quantum critical point. Upon 
reaching the quantum critical line from the above at 
T — > the system undergoes the first order quantum 
phase transition, which is FCQPT taking place at T c = 0. 
While at diminishing temperatures, the systems located 
before QCP do not undergo a phase transition and their 
behavior transits from NFL to LFL. 

It is seen from Fig. [10] that at finite temperatures 
there is no boundary (or phase transition) between the 
states of systems located before or behind QCP shown 
by the arrow. Therefore, at elevated temperatures the 
properties of systems with x/x-pc < 1 or with x/x-pc > 1 
become indistinguishable. On the other hand, at T > 
the NFL state above the critical line and in the vicinity 
of QCP is strongly degenerate, therefore the degeneracy 
stimulates the emergence of different phase transitions 
lifting it and the NFL state can be captured by the 
other states such as superconducting (for ex amp le, by 
the superconducting state (SC) in CeCoIns p54ll70l |) or by 
antiferromagnetic (AF) state (e.g. AF one in YbRh2Si2 
[4l|) etc. The diversity of phase transitions occurring 
at low temperatures is one of the most spectacular 
features of the physics of many HF metals. Within the 
scenario of ordinary quantum phase transitions, it is 
hard to understand why these transitions are so different 
from one another and their critical temperatures are so 
extremely small. However, such diversity is endemic to 
systems with a FC [23j. 

As mentioned above, the system located above QCL 
exhibits the NFL behavior down to lowest temperatures 
unless it is captured by a phase transition. The behavior 
exhibiting by the system located above QCL is in 
accordance with the experimental observations that 
study the evolution of QCP in YbRh2Si2 under negative 
chemical pressure induced with Ir/Ge substitution [7ll 
[72|| . As reported, the negative pressure leaves an 
intermediate spin-liquid-type ground state over a finite 
magnetic field range. We assume a simple model that 
the application of negative pressure reduces x and the 
system moves to a position over QCL. That is the 
Ir/Ge substitution drives YbRli2Si2 to a paramagnetic 
state over the finite magnetic field B range, so that 
the electronic system of the HF metal is located above 
QCL shown in Fig. [TU] and demonstrates the NFL 
behavior. In that case magnetic field B represents 
an additional control parameter driving YbRh2Si2 to 
the paramagnetic state. We predict that at lowering 
temperatures, the electronic system of YbRli2Si2 is 
captured by a phase transition, since the NFL state 
above QCL is strongly degenerate. At diminishing 
temperatures, this degeneracy is to be removed by some 
phase transition which likely can be detected by the LFL 
state accompanying it. 



Upon using nonthermal tuning parameters like the 
number density x, the NFL behavior can be destroyed 
and the LFL one will be restored. In our simple model, 
the application of positive pressure P makes x grow 
removing the QCP of FCQPT shown in Fig. [TU1 so 
that the electronic system of YbRli2Si2 moves into the 
shadow area characterized by the LFL behavior at low 
temperatures. The new location of the system is shown 
by the arrow pointing at the solid square. As a result, 
the application of magnetic field B ~ B c q does not 
drive the system to its QCP with the divergent effective 
mass because the QCP is already destroyed by the 
pressure. Here B c q is the critical magnetic field that 
eliminates the corresponding AF order. At B > B c a 
and raising temperatures, the system moving along the 
vertical arrow transits from the LFL regime to the NFL 
one. As seen from Fig. [TOl in the NFL regime the system 
properties coincide with that of the system at P = 0. 
The observed behavior is in accord with the experimental 
facts indicating the separation of the AF order from QCP 
under the application of pressure in YbRh2Si2 (73| . 



6. HIGH MAGNETIC FIELDS 
THERMODYNAMICS OF HEAVY FERMION 
QUASIPARTICLES 

As was mentioned in Section ® one of the most 
interesting and puzzling issues in the research of HF 
metals is the NFL behavior taking place in wide range 
of both temperature T and magnetic field B variations 
[73 |. For example, measurements of the specific heat on 
YbRh.2Si2 under the application of magnetic field B show 
that the range extends at least up to twenty Kelvins 
as shown in Fig. [3J In Fig. [5] the normalized effective 
mass as a function of normalized temperature T/y 
is shown by geometrical figures. Figures [3] and [6] reveal 
the remarkably large temperature ranges over which the 
NFL behavior and scaling one are observed. 

In this Section, we analyze the thermodynamic 
properties of YbRli2Si2 in both low and high magnetic 
fields. Our calculations of the specific heat and 
magnetization allow us to conclude that under the 
application of magnetic field the heavy-electron system of 
YbRli2Si2 evolves continuously without a metamagnetic 
transition. At low temperatures and high magnetic 
fields B ~ B* the system is completely polarized 
and demonstrates the LFL behavior, at elevated 
temperatures the HF behavior and related NFL one are 
restored. The obtained results are in good agreement 
with experimental facts as both temperature and 
magnetic field vary by more than two orders of 
magnitude. 

We use the Landau equation (|14l) to study the behavior 
of the effective mass M*(T,B) as a function of the 
temperature and the magnetic field. For the model 
of homogeneous HF liquid at finite temperatures and 
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magnetic fields, Eq. (jT^J) acquires the form 
1 1 f PfP 



M*(T,B) 



E 



F 

3 x cr,(Ti 



dn ai (p,T,B) dp 



>F,P) 



9p 



(2tt) 3 



(26) 



where F a ^ ai (pp, p) is the Landau amplitude dependent 
on momenta p-p, p and spin a. For the sake of defmiteness, 
we assume that the HF liquid is 3D liquid. The 
quasiparticle distribution function can be expressed as 



n CT (p,T) = I 1 + cxp 



(e CT (p,T) - fi a 



T 



(27) 



where e cr (p,T) is the single-particle spectrum. In our 
case, the chemical potential may depend on the spin 
due to the Zeeman splitting fj, a = fj, ± (J-bB, where 
I^ib is the Bohr magneton and /i is the chemical 
potential preserving the number of particles. We write 
the quasiparticle distribution function as n <T (p,T, B) = 
Sn a (p,T,B) +n CT (p,T = 0,B = 0). Equation dHJ) then 
becomes 



1 



M*(T,B) 
F a , ai (pf,Pi) 
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w* 



PfPl 

Pf 



d<5n CT1 (pi,T,£) dpi 



9pi 



(2tt) 3 



(28) 



Equation (|28p is solved with a quite general form 
of Landau interaction amplitude [75]. Choice of the 
amplitude is dictated by the fact that the system has to 
be at the quantum critical point of FCQPT, which means 
that the first two p-derivatives of the single-particle 
spectrum e(p) should equal zero. Since the first derivative 
is proportional to the reciprocal quasiparticle effective 
mass 1/M*, its zero just signifies QCP of FCQPT. The 
second derivative must vanish; otherwise e(p) — [i has the 
same sign below and above the Fermi surface, and the 
Landau state becomes unstable [H, [23j. Zeros of these 
two subsequent derivatives mean that the spectrum e(p) 
has an inflection point at pp so that the lowest term of 
its Taylor expansion is proportional to (p — pf) 3 - 

When the system is near FCQPT, it turns out that 
the normalized solution of Eq. (f2"8"]) M^(y) can be 
well approximated by a simple universal interpolating 
function. The interpolation occurs between the LFL 
(M* (xa + bT 2 ) and NFL (M* cx T~ 2 / 3 ) regimes!! 



M* N {y) « c 



l + c lV 2 



1 



c 2 y 



8/3 ' 



(29) 



Here a and b are constants, cq = (1+C2)/(1+Ci), c\ and C2 
are fitting parameters, related to the Landau amplitude. 
Magnetic field B enters Eq. (|27|) as B^lb/T making Tm oc 
B[ib [1, IzHl- We conclude that under the application of 
magnetic field the variable 



y = T /Tm oc 



T 



[i B {B-B c 



(30) 



remains the same and the normalized effective mass is 
again governed by Eq. (|29|) . Here B c q is the critical 
magnetic field driving both HF metal to its magnetic 
field tuned QCP and corresponding Neel temperature to 
T = 0. In some cases as in the HF metal CeRu2Si2, 
B c n = 0, see e.g. [zl|, while in YbRh 2 Si 2 , B c0 ~ 0.06 
T [27j and in CeCoIns the critical field is much larger 
B c q ~ 5.0 T [3l[. Therefore, in our simple model of the 
HF liquid B c q is taken as a parameter. In what follows, we 
compute the effective mass using Eq. (|26[) and employ Eq. 
for qualitative estimations of the considered values. 
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Phc. 11: Schematic phase diagram of YbRti2Si2 based on 
Eq. (|29|) with magnetic field as the control parameter B. 
The vertical and horizontal arrows show LFL-NFL and NFL- 
LFL transitions at fixed B and fixed T correspondingly. At 
B < Bco the system is in antiferromagnetic (AFM) state. 
The width of the transition regime T w oc T is shown by 
the segment confined by the two vertical arrows. Inset shows 
a schematic plot of the normalized effective mass versus 
the normalized temperature. Transition regime, where M% 
reaches its maximum value at y = T/Tm = 1, is shown by 
the hatched area both in the main panel and in the inset. 

We now are in position to construct the schematic 
phase diagram of the HF metal YbRh 2 Si 2 at B B* . 
The pase diagram is reported in Fig. [TT] The magnetic 
field B plays a role of the control parameter, driving 
the system towards its QCP. In our case this QCP 
is of FCQPT type. The FCQPT peculiarity occurs at 
B = B c o, yielding new strongly degenerate state at 
B < B c q. To lift this degeneracy, the system forms 
either superconducting (SC) or magnetically ordered 
(ferromagnetic (FM), antiferromagnetic (AFM) etc) 
state s |8| . In the case of YbRli2Si2, this state is AFM 
one [38]]. As it follows from Eqs. (|2"9"|) and (|50|) and 
seen from Fig. [TTJ at B > B c q the system is in 
either NFL or LFL states. At fixed temperatures the 
increase of B drives the system along the horizontal 
arrow from NFL state to LFL one. On the contrary, at 
fixed magnetic field and raising temperatures the system 
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transits along the vertical arrow from LFL state to NFL 
one. The inset to Fig. [TT] demonstrates the behavior of 
the normalized effective mass = M* /M M versus 

normalized temperature y = T/Tm following from 
Eq. (HU). The T~ 2 / 3 regime is marked as NFL one 
since (contrary to LFL case where the effective mass 
is constant) the effective mass depends strongly on 
temperature. It is seen that temperature region y ~ 1 
signifies a transition regime between the LFL behavior 
with almost constant effective mass and NFL one, given 
by T~ 2 / 3 dependence. Thus, temperatures T ~ Tm, 
shown by arrows in the inset and main panel, can be 
regarded as a transition regime between LFL and NFL 
states. It is seen from Eq. (|30f that the width of the 
transition regime T w oc T is proportional to (B — B c q). It 
is shown by the segment between two vertical arrows in 
Fig [TTJ These theoretical results are in good agreement 
with the experimental facts [H, [77} • 

In Fig. our calculations of the normalized effective 
mass A4^(Tn) at fixed magnetic field B* making the 
system fully polarized are shown by the solid line. 
Figure [T^] reveals the scaling behavior of the normalized 
experimental curves - the scaled curves at different 
magnetic fields B merge into a single one in terms of 
the normalized variable y = T /Tm- It is seen from Fig. 
1121 that our calculations are in accord with facts (H|: At 
elevated temperatures (y ~ 1) the LFL state first changes 
into the transition regime and then disrupts into the NFL 
state. 
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Phc. 12: The normalized effective mass — M* /M% f versus 
normalized temperature Tjv = T/Tm. is extracted from 
the data shown in Fig. [3] Magnetic fields B is listed in 
the legend. Our calculations (made at B ~ B* when the 
quasiparticle band is fully polarized) are depicted by the solid 
curve tracing the scaling behavior of Mjy. 



To get further insight into the behavior of the system 
at high magnetic fields, in Fig. [T31 we collect the curves 
Af^(Tjv) both at low (symbols in the upper box in Fig. 



ITS")) and high (symbols in the lower box) magnetic fields 
B. All curves have been extracted from the experimental 
facts [13, Hl| . It is seen that while at low fields the low- 
temperature ends (TV ~ 0.1) of the curves completely 
merge, at high fields this is not the case. Moreover, 
the low-temperature asymptotic value of C/T = M% 
at low fields is around two times more then that at 
high fields. The physical reason for low-field curves 
merging is that the effective mass does not depend 
on spin variable so that the polarizations of subbands 
with cr-j- and aj. are almost equal to each other. This 
is reflected in our calculations, based on Eq. (|2"5)) fol- 
low magnetic fields B « B* . The result is shown 
by the dotted line in Fig. [T3] It is also seen from 
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Phc. 13: Joint behavior of the normalized effective mass 
M^f at low (upper box symbols) and high (lower box 
symbols) magnetic fields extracted from the specific heat 
(C/T) measurements of the YbRh 2 Si 2 [H- Our low-field 
calculations are depicted by the dotted line tracing the scaling 
behavior of M^j. Our high-field calculations (solid line) are 
taken at B ~ B* when the quasiparticle band becomes fully 
polarized. 

Fig. Q2] that all low-temperature differences between 
high- and low field behavior of the normalized effective 
mass disappear at high temperatures. In other words, 
while at low temperatures the values of for low 
fields are two times more then those for high fields, at 
temperatures Tjf > 1 this difference disappear. It is seen 
that these high temperatures lie about the transition 
region, marked by hatched area in the inset to Fig. [TT] 
This means that two states (LFL and NFL) separated by 
the transition region are clearly seen in Fig . [TBI displaying 
good agreement between our calculations (dotted line 
for low fields and thick line at high fields) and the 
experimental points (symbols). 

Contrary to the low fields case, at high fields B ~ 
B* (symbols in the lower box of Fig. [TBJ, in the low 
temperature LFL state the curves M^(Tn) do not merge 
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into single one and their values decrease as B grows 
representing the full spin polarization of the HF band 
at the highest reached magnetic fields. As we have 
mentioned above, at temperature raising all effects of 
spin polarization smear down, yielding the restoration 
NFL behavior at T ~ [ibB. Our high-field calculations 
(solid line in Fig. [T3"]) reflect the latter fact and are also 
in good agreement with experimental facts. In order not 
to overload Fig. [T3] with unnecessary details, we show the 
calculations only for the case of the full spin polarization. 
Figure [HI reports the maxima M^(B) of the functions 
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Phc. 14: The maxima M M (B) of the functions C/T versus 
magnetic field B for YbRli2Si2. The maxima taken at the 
low fields [37l ] are shown by the solid squares, the ones taken 
at the high fields [38j are depicted by the solid circles. The 
solid curve is approximated by M M (B) oc d/s/~B~- 
a fitting parameter. 



Bdo, d is 



C/T shown in Fig. Inversus B depicted by the solid circles 
and the solid squares, corresponding measurements on 
YbRli2Si2 taken at the low magnetic fields [37]. The 
solid line represents the approximation describing the 
maxima value of M*(B) oc \JB — B c n calculated within 
the framework of FCQPT theory % It is seen 

that our calculations are in good agreement with the 
experimental facts over two decades in magnetic field B. 
The agreement indicates that at T ~ [ibB the transition 
regime takes place and the NFL behavior restores at 
higher temperatures T ~ 20 K. 

In Fig. Qj)l the solid squares and circles denote 
temperatures Tn(B) at which the maxima of C/T take 
place. To fit the experimental data |37t IHJ the function 
T N (B) = b(B - B c0 ) defined by Eq. ||30J| is used. It is 
seen from Fig. [15] that our calculations shown by the 
solid line are in accord with experimental facts, and 
we conclude that the transition regime of YbRh 2 Si2 is 
restored at temperatures T ~ [IbB. Consider now the 
magnetization M(B, T) as a function of magnetic field B 
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Phc. 15: The temperatures Tm(B) at which the maxima of 
C/T are located versus magnetic field B for YbRli2Si2. The 
solid squares denote Tm taken in measurements at the low 
magnetic fields [3711 - The solid circles denote Tm at the high 
magnetic fields [3a |. The solid line represents the function 
Tm oc b(B — B c q), b is a fitting parameter, see Eq. (|30[> . 



at fixed temperature T = Tt 



M(B,T)= / X (z,T)dz 



(31) 



where the magnetic susceptibility x is given by [l5l [l6| 

j3M*(B,T) 



X(B,T) 



1 



(32) 



Here, ft is a constant and Fft is the Landau amplitude 
related to the exchange interaction. 

At low magnetic fields, B <C B* , our calculations 
show that the magnetization exhibits a kink at some 
magnetic field B — Bk- The experimental magnetization 
demonstrates the same behavior [3^, [7{|. We use St- 
and M(Bfc) to normalize B and M respectively. Due to 
the normalization the coefficients /3 and (1 + Fft) drop 
out from the result, and x <x M* HI- So we can use 
Eq. (|28p to calculate the magnetic susceptibility x- The 
normalized magnetization M(B)/M(Bk) both extracted 
from experiment (symbols) and calculated one (solid 
line), are reported in the inset to Fig. [TCJ It shows 
that our calculations are in good agreement with the 
experimental facts, and all the data exhibit the kink 
(shown by the arrow) at Bjy ~ 1 taking place as soon as 
the system enters the transition region corresponding to 
the magnetic fields where the horizontal arrow in Fig. [TT] 
crosses the hatched area. To illuminate the kink position, 
in the Fig. [16] we present the M(B) dependence in 
logarithmic - logarithmic scale. In that case the straight 
lines show clearly the change of the power dependence 
M{B) at the kink. At magnetic field B ~ B* the 



15 




Normalized field 



1E-5 ' " lE-4 ' "lE-3 0,01 
Magnetic field B, arbitrary units 



Phc. 16: The calculated magnetization M(B) (symbols) and 
straight lines which are guides for eye. The intersection of the 
straight lines visualize the kink at the crossover region in the 
Fig. 1111 The inset: The field dependencies of the normalized 
magnetization M of YbRu 2 Si 2 [H at T = 0.05 K. The kink 
(shown by the arrow) is clearly seen at the normalized field 
Bn = B/Bk — 1. The solid curve represents our calculations. 



quasiparticle band becomes fully polarized and a new 
kink appears (38l [80| . we call this kink as the second 
one. Our calculations of the normalized magnetization 
(line) and the experimental points (squares) are shown 
in Fig. [T7] In that case both the magnetization and 
the field are normalized by the corresponding values at 
the second kink. In the Fig. [T7] we plot our calculated 
magnetization (in the variables B and M normalized 
to those in the second kink point) versus experimental 
one. The marvelous coincidence is seen everywhere except 
the high-field part at Bn > 1 where the normalized 
magnetization M nor exhibits a linear dependence on 
Bn depicted by the two arrows, while the calculated 
magnetization is approximately constant at these Bn- 
Such a behavior is the intrinsic shortcoming of the HF 
liquid model that takes into account only heavy electrons 
and omit the other conduction electrons that contribute 
to the magnetization (8ll. |82|. Thus, we can consider the 
high- field (at Bn > 1) part of the magnetization as 
the contribution which does not included in our theory. 
To separate this contribution from the experimental 
magnetization curve, we (numerically) differentiate it, 
then subtract constant part at Bn > 1 and integrate 
back the resulting curve. The coincidence between our 
calculations depicted by the solid curve and processed 
experimental data shown by the stars is reported in the 
inset to Fig. [TTJ As we can see now, the coincidence 
between the theory and experiment is good in the entire 
magnetic field domain. Taking into account the obtained 
results displayed in Figs. [T5HT71 we conclude that the 




Phc. 17: The normalized magnetization M nor as a function 
of the normalized magnetic field Bn- The calculated 
magnetization is shown by the solid line, and the experimental 
facts [38| is represented by the solid squares. The linear 
dependence of M nor on Bn is shown by the arrows. The 
inset demonstrates the experimental data with the subtracted 
linear dependence depicted by the stars. Our calculations is 
shown by the solid line. 



HF system of YbRu2Si2 evolves continuously under the 
application of magnetic fields. This observation is in 
agreement with experimental observations [Hj|. 



7. SUMMARY 

We have shown that the physics of systems with heavy 
fermions is determined by the extended quasiparticle 
paradigm. In contrast to the paradigm that the 
quasiparticle effective mass is a constant, within the 
extended quasiparticle paradigm the effective mass of 
new quasiparticles strongly depends on the temperature, 
magnetic field, pressure, and other parameters. The 
quasiparticles are well defined and can be used to describe 
the scaling behavior of the thermodynamic HF metals in 
a wide range of both temperatures T and magnetic fields 
B. 

It was demonstrated that the phase diagram of systems 
located near the fermion condensation quantum phase 
transition are strongly influenced by control parameters 
such as chemical pressure, pressure, the number density 
or magnetic field. Analyzing the phase diagram of 
YbRli2Si2, we have found that under the application 
of the chemical pressure (positive/negative) QCP is 
destroyed or converted into the quantum critical line, 
correspondingly. 

We have analyzed the thermodynamic properties of 
YbRli2Si2 in both the low and the high magnetic 
fields. Our calculations allow us to conclude that in 
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magnetic field the HF system of YbRli2Si2 evolves 
continuously without a metamagnetic transition and 
possible localization of heavy 4/ electrons. At low 
temperatures and under the application of magnetic 
field the HF system demonstrates the LFL behavior, at 
elevated temperatures the system enters the transition 
region followed by the NFL behavior. 

Our observations are in good agreement with 
experimental facts and show that the fermion 
condensation quantum phase transition is indeed 
responsible for the observed NFL behavior and 



quasiparticles survive both high temperatures and 
high magnetic fields. 
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